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ABSTRACT 

The coronae in X-ray binaries (XRBs) still remain poorly understood, although they have been 
believed for a long time to play a key role in modeling the characteristic outbursts of XRBs. Type-I 
X-ray bursts, the thermonuclear flashes happening on the surface of a neutron star (NS), can be used 
as a probe to the innermost region of a NS XRB, where the corona is believed to be located very 
close to the NS. We report the discovery of a tiny life cycle of the corona that is promptly co-evolved 
with the type-I bursts superimposed on the outburst of the NS XRB IGR J17473— 2721. This finding 
may serve as the first evidence of directly seeing the rapid disappearance and formation of a corona 
in an XRB with a cooling/heating timescale of less than a second, which can strongly constrain the 
accretion models in XRBs at work. 

Subject headings: stars: coronae — stars: neutron — X-rays: individual(IGR J17473-2721) — X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 

The spectra of low-mass X-ray binaries (LMXBs), in- 
cluding both black hole (BH) and neutron star (NS) 
XRBs, are generally disentangled into soft/thermal (e.g., 
a blackbody-like) and hard/Comptonized (e.g., power- 
law shape with a cutoff at tens to hundreds keV) com- 
ponents. For example, a hard X-ray power law with 
a cutoff energy of tens keV in a NS LMX B was first 
obser ved in the atoll source 4U 1608—522 (jZhang et al.l 
1996). A corona in a region close to the compact ob- 
ject of an XRB is usually invoked to account for the 
hard X-rays produced via inverse Compton scattering of 
the soft photons off the hot electrons in the corona. Al- 
though the existence for such a corona has been hinted 
from time to time, so far a direct observational evi- 
dence of its existence is still missing. It also remains 
a puzzle on how the corona is formed and where it is 
located. Plasmas may be energized to form a corona 



via e ither evap orati on dMever et al.1 [19941 : lEsin et ahl 
119971 iLiu et alJl2007t iFrank et al.H2002D or magnetic re- 



conn ections ([Zhang et alj|2000t lZhandl2007tlMever et al.l 
2007). The corona's location is sti ll not known for NS 
XRBs, and accordingly the Eastern (Mitsuda et al.lll989fl 
and the Western models ([White et all 1198 8) were pro- 
posed. The Eastern model considers a corona around 
the NS, whereas the Western model takes t he corona 
cover ing the accretion disk. A hybrid model ([Lin et al.l 
2007) for NS X-ray binaries was devised by analyzing 
the X-ray spectra of Aql X-l and 4U 1608-52, and of- 
fered a weak-Comptonization solution that differs from 
the above two models. X-ray binaries (XRBs) are com- 
monly observed with outbursts displaying strong spec- 
tral evolution. Here we take advantage of the type-I 
X-ray bursts from the hard surface of a NS to probe 
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the purported corona, under a circumstance of out- 
burst and accompanying spectral evolution. Type-I X- 
ray bursts are caused by unstable burning of the ac- 
creted hydrogen/helium on the surface of a NS, and 
manifest themselves as a sudden increase (typically 
by a factor of 10 or greater) in the X-ray luminos- 
ity (for reviews, see (iLewin et alj|1993t iCumming 120041 : 



IStrohmaver fc Bildstenll2006t iGallowav et alJ l2008) ) 



2. OBSERVATIONS AND DATA ANALYSIS 

IGR J 17473 -2721 was discover ed by INTEGRAL dur - 
ing an outburst in April 2005 (Gr ebenev et al.l 12005). 
Its subsequent outburst in 2008 March - October was 
almost uniformly monitored by RXTE/PCA (Propor- 
tional Counters Array), resulting in 182 RXTE/PCA 
pointed observations, with the identifier (OBSID) of pro- 
posal number (PN) 93064, 93093, and 93442 in the High 
Energy Astrophysics Science Archive Research Center 
(HEASARC). These observations adopted in this analy- 
sis cover the entire outburst and add up to ~ 475 ks of 
exposure time on the source. The analysis of the PCA 
data is performed by using HEAsoft v. 6.6. The data are 
filtered using the standard RXTE/PCA criteria. Only 
the data from the PCU2 (the third Proportional Coun- 
ters Unit, in the 0-4 numbering scheme) are used for the 
analysis, because this PCU was 100% working during all 
the observations. The background file used in the analy- 
sis of PCA data is the most recent one for bright-sources 
found at the HEASARC website 3 . The dead time cor- 
rection is made to all the spectra and lightcurves under 
the standard procedure described at the HEASARC web- 
site 4 . 

3. RESULTS 

As shown in Figure [1] the source in the 2008 out- 
burst experienced a two-months preceding low/hard 
state (LHS) and a lagging LHS with respect to the 
high/soft state (HSS); here 'low' and 'high' refers to 

3 pca_bkgd_cmbrightvle_eMv2005 1128. mdl 

4 http://heasarc.nasa.gov/docs/xte/recipes/pca_deadtime.html 
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its flux in the 2-10 keV band, and 'soft' and 'hard' 
refers to its spectral shape from 2-50 keV. The pre- 
ceding LHS and the lagging LHS are usually refereed 
in literature as the high lumi nosity LHS and the low 
lumin osity LHS, respectively (jRcmillar d fc McClintockl 
2006). The flux of the lagging LHS is roughly 
4 times lower than the preceding LHS, presenting 
with a so-c alled hysteresis typica l to the outbursts of 
BH XRBs (iMivamoto et all [1995t iMaccarone fc Coppil 



20031: iMever-H ofmciste r et all 120051 : iDunn et all l2010t 
Tang et al 11201 lh. Besides 4U 1908+ 005 (Aql X-l) and 
4U 1608-522~ i|Gladstone et all[2007h . IGR J17473-2721 
is thus the third NS XRB observed with hysteresis. 
There are 32 bursts located in the preceding LHS and 
10 bursts in the tail of the HSS and the lagg ing LHS 
(jChen et all [201(1 I20TT1 : IChenevez etaD l20ll . consti- 
tuting an entire sample of 42 bursts recorded by RXTE. 
The long-lived preceding LHS and the large number of 
burst events occurred in both the preceding and lagging 
LHSs make IG R J 17473 -2 721 a lmost unique for prob- 
ing its corona (I Zhang et al.ll2009() under a circumstance 
of an outburst evolution. 

For both LHSs, the persistent emission can be fitted 
by a power-law model with a cutoff energy ~40 keV. 
The persistent luminosity in the preceding LHS and the 
lagging LHS are ~0.1 £Edd (Eddington luminosity) and 
^0.025 Le<m, respectively. The burst spectrum is well 
modeled by a blackbody model of a characteristic tem- 
perature of less than 3 keV. The burst emission can reach 
^Edd> and dominates the total emission at energies well 
below ~30 keV, above which the persistent emission from 
the purported corona dominates. We therefore take the 
emission in the 30-50 keV energy band to investigate the 
possible influence of the bursts upon the corona. 

We choose 40 bursts out of the entire sample recorded 
by RXTE in the 2008 outburst with complete data for 
both the burst and the persistent emission. For each 
burst, we use the time when the burst reached its peak 
at 2-10 keV as a reference to produce the lightcurve and 
spectrum of each burst. Those X-rays recorded 48 sec- 
onds before and 80 seconds after the reference time are 
regarded as the background and are subtracted off for 
each burst in timing and spectral analysis. The 30-50 
keV persistent count rate recorded by RXTE/PCA in the 
preceding LHS, the lagging LHS and the HSS are ~4.0 
cts/s, ~1.0 cts/s, and ~0.11 cts/s, respectively. With 
respect to the reference time we split each burst into 
a sequence of 16-second segments. After the persistent 
emission is subtracted off, bursts are combined for those 
located in the preceding LHS or in the lagging LHS. As 
shown in Figure [2] the 30-50 keV flux of the combined 
burst in the preceding LHS is mostly negative during the 
burst and around zero elsewhere; however, the fluxes of 
the combined burst in the lagging LHS are always con- 
sistent with zero. 

Our previous reports showed that the corona of the 
system (Figure [3]) is most likely located around a largely 
inclined disk and forms a relatively small opening angle 
with respect to the surface of neutron star, so that the 
type-I X-ray bursts can escape from t he NS surface with- 
out suffering seve re Comptonization (|Zhang et all 120091 : 
iChen et al.ll201ll ). The 30-50 keV decrement reaches a 
maximum of ~2 cts/s at the 2-10 keV burst peak, ac- 
counting to about half of the 30-50 keV persistent flux. 



This suggests that half of the corona were cooled by the 
soft photons of the bursts. Actually the 30-50 keV profile 
in the preceding LHS is anti-correlated with that of the 
2-10 keV under a correlation coefficient of 0.89 (see Fig- 
ure [4|). In this analysis, a time resolution of 16 seconds 
is used to show the significance of the flux decrement, 
due to the very low flux in the 30-50 keV band. We try 
as well to extract the net lightcurves of the combined 
burst in the 2-10 keV and 30-50 keV bands with a time 
resolution of 1 second in the preceding LHS. A cross- 
correlation analysis between the two lightcurves shows 
that the 30-50 keV X-rays lag the 2-10 keV X-rays by 
0.7±0.5 second, as shown in Figure [SJ This suggests 
that the fading and recovering of the corona follow the 
burst flux change almost instantaneously. In the preced- 
ing LHS, the Compton cooling time of a pure electron- 
positron plasma is i com pt = E c /P compt = 5.4 x 10~ 8 s, 
here E c ~ 20 keV is the electron energy. However, the 
plasma is more likely a proton-electron plasma, whose 
cooling time is ~ 10 -7 s, i.e., increased by a factor of 
2. This means that the heating time of the corona must 
be longer than ~ 10 -7 s; otherwise the corona cannot 
be cooled down as observed during the rising phase of 
the bursts. On the other hand, the heating time must 
be shorter than about 1 second; otherwise the recovering 
of the corona would lag the bursts in the falling phase. 
These results serve as the first observational evidences 
in an XRB for a promptly cooling and recovering corona 
which manifests itself with an apparent tiny life cycle via 
tightly co-evolving with the soft X-ray bursts. 

It is normally difficult to estimate the mass of the 
corona in an XRB, because the seed photon energy den- 
sity C/ph for the inverse Compton scattering in the corona 
cannot be observed directly and cleanly. Fortunately, the 
prompt response of the corona to the bursts offers an 
unique opportunity to estimate the mass of the corona 
directly. The temperature and Compton scattering opti- 
cal depth of the corona are obtained to be around 20 keV 
and 2 for both LHSs, respectively, by fitting their persis- 
tent emission spectra with the thermal Com ptonization 
model (Comptt in Xspec) (jTitarchukl 11994} ) . With the 
bursts as the coolant, the seed photon energy density is 
give by, 

u * - iSt' (1) 

here R is the distance of the corona from the NS. The 
total power for a single scattering from a non-relativistic 
electron with temperature fcT, i.e., the cooling rate, is 
given by, 



AkT 



compt 



( j)ccr T £/p h (erg s 1 ), 



mc 



(2) 



here ctt and U p h are the Thomson cross section and pho- 
ton energy density of a burst in the corona, respectively. 
The mass of the corona cooled by a burst can be calcu- 
lated as, 

ALlhs ,„s 

-m H . (3) 



AM, 



corona 



P 



comp 



The inner radius of the accretion disk is assumed as the 
radius of the magnetosphere of the NS, where the gas 
pressure in the disk is balanced by the magnetic pres- 
sure of magnetosphere (Cui 1997; W eng fc Zhandl2011j) . 
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Taking L burst = 8.8 x 10 37 erg s l ~ 0.5 i E dd , ^lhs = 
1.8 xlO 37 ergs" 1 - 0.1 i E dd, the mass of NS - 1.4 M©, 
the NS surface magnetic field strength ~ 10 8 G, the 
radius of the NS ~ 10 km, we have, R — 12 km and 
Mcorona = 2AAf corona = 5.8 x 10 13 g. Since the Comp- 
ton scattering optical depth is around 2 and in the above 
calculation a single scattering is assumed, the calculated 
corona mass is only an upper limit. The accretion rate 
in the preceding LHS M = %P = 1.0 X 10 17 g s _1 , 
here e = 0.2 is the radiative efficiency for a NS XRB. 
Therefore the accretion process has sufficient mass reser- 
voir to supply the required mass to the corona almost 
instantaneously, as observed. 

4. DISCUSSION 

Despite the similar spectral properties of the persis- 
tent emission for both the preceding and lagging LHSs, 
no decrement of the 30-50 keV photons is observed in the 
combined burst in the lagging LHS. The lower flux in the 
lagging LHS indicates an accretion rate decreasing by a 
factor of 4 compared to that in the preceding LHS. There- 
fore the inner disk radius (R oc L~ 2//7 ) (|Weng fc Z hang 
2011) in the LHS is almost doubled, resulting a lower 
cooling rate by the bursts. With the persistent flux of 
about 1 cts/s in the lagging LHS in the 30-50 keV band, 
the maximum of decrement in the combined burst is ex- 
pected to be ~0.2 cts/s, well within the observed uncer- 
tainties. This explains the lack of flux decrement in the 
30-50 keV band in the lagging LHS. 

Spectral transitions in XRB outbursts have been stud- 
ied in the frameworks of producing/supressing the hard 
X-rays via mechanisms, e.g., disk viscous dissipation and 
Compton cooling, where one key parameter is the time 
scale of flux variability. While the inverse Compton scat- 
tering has been applied successfully to account for the 
cooling, however no consensus has been reached on the 



correct corona formation model, e.g., between disk evap- 
oration and magnetic re-connection models. Unfortu- 
nately, the time scales for either the transition into a 
HSS or stepping back to a LHS have so far only been 
observed not shorter than days, not sufficiently fine to 
discriminate between different models. For the disk evap- 
oration, the formation of a corona is driven and energized 
by th e disk accretion, which has a typical time scale of 
davs rtMever et aLl ll994t lEsin et al.lll997t iLiu et al.ll2007t 
iFrank et al.l l2002t ) , obviously inconsistent with the ob- 
served heating time of shorter than 1 second. On the 
other hand, magnetic field reconnections in the inner disk 
region can release the kinetic energies in the rotating disk 
to heat the corona with the Keplerian orbital time scale 
of the order of milliseconds. 

Therefore microscopic processes in the accretion disk 
should be responsible for forming the coronae in this NS 
XRB, in a process similar to that powering the solar 
flares in the Sun, a s already discussed previously for BH 
accreting sy stems (|Zhang et al.ll2000t IMever et all 120071; 
Zhang 2007). Future X-ray instruments with fine tim- 
ing resolution and much larger effective areas than the 
RXTE/PCA instrument, may allow direct probing the 
fine details of the microscopic processes in the accretion 
and corona formation, by observing the interplays be- 
tween type I X-ray bursts and the rapid spectral evolu- 
tion of NS XRBs. ' 
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Fig. 1 — RXTE/PCA lightcurve (2-10 keV, upper panel) and Swift/BAT lightcurve (15-50 keV, lower panel) covering the 2008 outburst 
of IGR J17473— 2721 with a time resolution of 1 day. The dashed lines show the different states of the outburst. The locations of the bursts 
are marked at the top of the upper panel (blue vertical lines). 
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Fig. 2. — The 16s-bin lightcurves for the bursts in the preceding LHS (upper panel) and the lagging LHS (lower panel). Each data point 
is the sum over the spectral residual after subtracting off the persistent emissions at 2-10 keV (black) and 30-50 keV (red), respectively. 
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Fig. 3. — The illustration of the central region of a NS XRB, in which a corona is located around the disk and cooled by the soft X-rays 
from a type-I burst occurred on surface of the NS. 
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Fig. 4. — A linear fit on the data shown in the upper panel of Figure 2. 
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Fig. 5. — The cross-correlation between the 2-10 keV and 30-50 keV, with a time resolution of 1 second, for the combined burst in the 
preceding LHS. 



